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We consider a growing network in which an incoming node gets attached to the j"* existing node 
with the probabihty Hi oc k^r", where ki is the degree of the i' node and n its present age. The 
phase diagram in the q — /3 plane is obtained. The network shows scale-free behaviour, i.e., the 
degree distribution P(k) ~ fc~^ with 7 = 3 only along a line in this plane. Small world property, on 
the other hand, exists over a large region in the phase diagram. 
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Complex web-like structures describe a wide variety 
of systems of high technological and intellectual impor- 
tance. The statistical properties of many such networks 
have been studied recently with much interests. Such net- 
works, with complex topology are common in nature and 
examples include the world wide web, the Internet struc- 
ture, social networks, communication networks, neural 
networks to name a few [1-3] . 

Some of the common features which are of importance 
in these networks of diverse nature are: 

(i) Diameter of the network: This is defined as the 
maximal shortest path over all vertex pairs. The net- 
works in which the diameter (P) varies as the logarithm 
of the number of nodes (N), i.e, V oc ln{N), are said to 
have the small world (SW) property. On the other hand 
if T) varies as a power of N we get what can be termed 
as regular networks. One can also study D, the average 
shortest distance between pairs of nodes. In general D 
and V have the same scaling behaviour. 

(ii) Clustering coefficient: A common property of many 
real networks is the tendency to form clusters or trian- 
gles, quantified by the clustering coefficient. This is 
in contrast to random networks (Erdos and Renyi [4]) 
where small world property is present but the clustering 
tendency is absent. 

(iii) Degree Distribution: The node degree distribu- 
tion function P{k) gives the probability that a randomly 
selected node has exactly k edges. In a random network 
this is a Poisson distribution. In many real world net- 
works -P(fc) shows a power law decay and such networks 
are termed Scale Free Networks. 

In order to emulate the different features of real 
networks several models have been proposed. While 
the Watts-Strogatz [5] network provides an appropriate 
model for the small world network (i.e., small diameter 
and finite clustering coefficient), scale free properties of 
a network can be reproduced by models proposed later 
by Barabasi and Albert (BA) [6] and independently by 



Hubermann and Adamic [7]. 

In the BA model, a network is grown from a few nodes 
and new nodes are added one by one. At a time i, the 
incoming node is connected randomly to the «*'' existing 
node with the attachment probability 11^ (i) given by 



n.(t) - fc,(t). 



(1) 



where ki is the degree of the i*'' node. The degree distri- 
bution in the BA model shows the scale-free behaviour 



P{k) ~ k"^, 



(2) 



with 7 = 3. 

Following its introduction, several modifications in the 
BA model have been studied. A few of them are worth 
mentioning here in the context of the present paper. Non- 
linear dependence of the attachment probability on k, 
in the model designed by Krapivsky et al (KRL) [8], 
shows that the scale-free property exists only for the lin- 
ear dependence. This nonlinear model has been studied 
in much detail very recently in [9]. On the other hand, 
BA model on an Euclidean network [10,11] has also been 
considered in which the attachment probability has been 
modified as follows: 



n,(t) ~ h{tfi' 



(3) 



where I is the Euclidean length between the new and old 
nodes. A phase diagram in the (3 — 8 plane was obtained 
along with other interesting features. 

Another important modification in the BA model has 
been made by incorporating time dependence in the net- 
work [12,13]. In real life networks, a time factor may also 
modulate the attachment probability. In most of the real 
world networks, aging of the nodes usually takes place, 
e.g., one rarely cites old papers, or in social networks, 
people of the same age are more likely to be linked. Doro- 
govtsev and Mendes (DM) [12] studied the case when the 
connection probability of the new site with an old one is 
not only proportional to the degree k but also to the 
power of its present age, r, such that 



n,(i)~fc,(i)rf 



(4) 



and they showed both numerically and analytically that 
the scale free nature disappears when a < — 1. Here a 
governs the dependence of the attachment probability on 
the "age difference" of two nodes, i.e, for negative values 
of a, a new node will tend to attach itself to the younger 
nodes. Therefore for the extreme case a — > — oo, a new 
node will attach itself to its immediate predecessor while 
for the case a ^ oo, the oldest and a few very old nodes 
will get more edges. The time dependence presents a 
competing effect when a < but for a > 0, the older 
nodes get even more rich, similar to the rich gets richer 
effect. 

Encouraged by the existence of the various phase tran- 
sitions observed in the modified BA models, we have fur- 
ther generalised the time dependent BA network. Here 
we generate a network such that the attachment proba- 
bility is given by 



D 



H,(t) ^ h{t) 



Pn 



(5) 



We expect here that /3 7^ 1 will change the behaviour of 
the DM model as in [8]. The competing effect of a is 
able to destroy the scale free nature of the DM model 
(/3 — 1). The effect of a positive /3(> 1) and negative a 
could re- instate the scale free behaviour as in [11] and it 
is also possible to obtain a phase diagram in the a — j3 
plane. Formally eq (5) is analogous to eq (3). However 
here the nodes are chronological, i.e, the age of the ini- 
tial node is t at time t, that of the second node i — 1 and 
so on. In the Euclidean network on the other hand the 
coordinates of the nodes are uncorrelated. Moreover, the 
dimensionality plays an important factor in it. 
The known limits of this model are therefore 

1. /3 = 1, a = - BA network 

2. 13 any value, a = - KRL network 

3. P = l,a any value - DM network 

When a and /3 are both zero, we get a random growing 
network which shows an exponential decay of P{k) [14]. 

The network is generated in the usual manner where 
we start with a single node and at every time step the 
new node gets attached to one of the existing nodes with 
an attachment probability given by equation (5). 

We have considered nodes with a single incoming link 
such that there are no loops and the clustering coefficient 
is zero. Thus we focus our attention on the degree distri- 
bution and the average distance to study the small world 
and scale free behaviour. 

From eq (5) we predict that for any /3 as a — > -|-cx) a gel 
formation is expected when majority of the nodes tend 
to get attached to a particular node. On the other hand 
when a — > — cxD we expect a regular chain formation (in 
the time space) when each node gets attached to its im- 
mediate predecessor. The average shortest distance (D) 
for both a -^ -}-cxd and a — ^ — cxd are easy to calculate. 
When a —> —00, I? is given by 



Sf (fc(fc ^ 1) + (N ^ k){N ^ k + 1)) 
2N{N-1) ' 



^ (iV + l)/3. 

On the other hand, for large values of a, D has a finite 
value ~ 0(1). Hence it is natural to expect a transition 
from a small world behaviour to a regular behaviour as a 
is varied. In fact for [3 = 0, one can locate approximately 
the transition point using some simple arguments. 
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FIG. 1. The variation of the average shortest distance D 
with A'^ for various values of a at (a)/3 = 2 and (b)/3 = 0.5. 
In (a) the exponent A changes value from 1 to sharply as 
we go from top to bottom of the figure, corresponding to 
a = — 10, — 5, — 3, — 2, 0. In (b) A changes from 1 to a rela- 
tively small value as a = —10, —5, —2, from top to bottom. 

In analogy with [15], one can define here an "age dif- 
ference factor" Atjj = \ti — tj\ if the i*'' node of age ti 
and j*^ node of age tj are connected. If the network has 
been evolved till a time t{> 2), then for the incoming 
node we can write 



(At), 



i: 



'dT 



\> 



dT 



(6) 



For the small world property, the behaviour of (At) 
for large t should be studied. From cq(6), for large t, 
At r^ 0(1) for a < —2 and therefore there can be no 
small world behaviour for a < —2 for large networks. 
On the other hand for a > — 1, (At) ^ 0(i), from which 
one can expect SW property for a > — 1. We in fact find 
a small world to regular network transition at a = — 1 
numerically. 



The simulations have been made using a maximum of 
2000 nodes for studying small world properties and 4000 
nodes for determining degree distribution, using a maxi- 
mum of 1000 configurations. In the analysis of the small 
world characteristics, wc calculate D for the networks 
for different values of /3 and a. The D versus N curve 
is generally of the form D ^ N^, where the exponent A 
depends on a (see Fig. 1). 

In order to locate the transition to the small world 
(where A is either zero or has a very small value) we note 
the variation of A with a for different values of (3. We 
observe that for all values of (3, there is a sharp fall in A 
from unity to a very small value indicating a transition 
from regular to small world behaviour. The transition 
point shifts to a more negative value of a as we move 
from smaller to larger values of (3. Typical A vs. a plots 
are shown in Fig. 2. 
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FIG. 2. The variation of the exponent A with a for differ- 
ent values of (3 (from right to left of the figure, /3 = 0, 2, 3, 4). 
There is a sharp transition in the value of a from -|-1 to as 
network behaviour changes from regular to small world. 

Next we study the degree distribution P{k) for the 
network for several values of a and (3. For the regular 
chain limit (a — > — cxd), most of the vertices have degree 
2, while for the gel phase (a -^ +cxd), there will be a 
very high maximal degree and many leaves (i.e, nodes 
with degree=l). Thus the different phases can be iden- 
tified from the behaviour of P{k). First let us discuss 
the known case for (3 — 1. We find that P{k) has an 
exponential decay at a = —1 as in [12] and has scale-free 
(SF) behaviour for a > —0.5. The latter value does not 
agree with [12] and the possible reasons of discrepancy 
are discussed later. For other values of a, P{k) has a 
stretched exponential (SE) behaviour, i.e.. 



P{k) - exp(-aa;^), 



(7) 



where b depends on a. Allowing /3 to assume values 
greater than unity, we find that SF behaviour exists only 
for a specific value of a == a*, e.g. at /3 = 3 we obtain 
this behaviour at a = a* = 2.5 (Fig. 3). 
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FIG. 3. (a) The P{k) vs k plot in log-log scale for /3 = 3. 
Here, at a = q* = 2.5 (dashed line), there is scale free be- 
haviour, while for other values (e.g. a = —2.8, —2.3) this 
behaviour is lost. (b)P(k) vs k plot in log-linear scale for 
l3 — 0.8. Here exponential behaviour is observed at a = —0.9 
(dashed line), while for other values (e.g. a — —0.5,-2.0) 
stretched exponential nature is observed. 

For a > a*, we get a gel-like behaviour, while for 
a < a* , we again get a stretched exponential behaviour. 
The scale-free behaviour for /3 > 1 always occurs with 
7 = 3 as in the BA network. For /3 < 1, SF behaviour is 
not observed for any value of a. Here P{k) shows a SE 
behaviour as in eq (7). 

In fig. 4, we have shown the phase diagram in the a — /3 
plane. We have plotted the phase boundary between 
the SW and the regular network regions, the curve along 
which scale- free behaviour exists and the line along which 
6=1. The b = 1 line is not a phase boundary, but it 
has the interesting property that it has the behaviour of 
a random growing network albeit with non-zero values of 
a and (3. For negative values of a, aging can be regarded 
as a competitive phenomenon to preferential attachment 
to the extent that one recovers the random behaviour 
even for high values of (3 along this line. 

A small world network is not necessarily scale-free but 
a scale- free network is usually a small world barring some 
exceptional or artificial cases (e.g, if one considers N 
number of BA networks in a series, it is a scale- free but 
not a small world network). To the right of the scale- 
free line and above it, the network shows a gel- like be- 
haviour. Both the scale-free and gel regions do have the 
small world property, as expected, but there are finite re- 
gions in the phase diagram where the degree distribution 
is not a power law but of different types (e.g., exponential 
or stretched exponential) with small world behaviour. 
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FIG. 4. The phase diagram for the given network in the 
a-P plane. The small world (SW) regions with gel-like as well 
as stretched exponential (SE) behaviour, the regular chain re- 
gion and the scale free region are indicated. The network is 
scale free along the thinner solid line while the broken line is 
the phase boundary for SW-regular transition and these two 
lines merge along the thicker solid line. The dotted line is the 
one along which 6=1, i.e, where the network is random in 
nature. (All lines are guides to the eye.) 

In summary, we have generalised the BA network to 
include time-dependent or aging effects in the attach- 
ment probabihty (eq. (5)) such that both the time de- 
pendence and the degree dependence can be parametri- 
cally tuned. A phase diagram is obtained in the a — P 
plane, where a, (3 are the parameters governing the two 
factors respectively. We claim that this is the most gen- 
eralised network where both time dependence and degree 
dependence are incorporated in the preferential attach- 
ment. There is a quantitative disagreement in the tran- 
sition point a,t f3 — I as compared to [12] which may 
be because of the finite sizes considered here. The time 
and effort required to locate phase transition points are 
considerable and a finite size analysis has not been at- 
tempted therefore. Other results known earlier, e.g., gel 
formation beyond /3 > 1 for a — 0, exponential decay of 
P{k) for both a,f3 = etc. have been recovered in our 
simulations. Similar to the Euclidean network [10,11], 
the scale free behaviour is found to exist along a single 
line here. In fact, as regards the scale-free boundary, 
the present phase diagram is very much similar to that 
obtained in [11]. However, here the network belongs to 
the BA universality class (7 = 3) along the entire line. 
Moreover, one can compare the present results with the 
one dimensional Euclidean network only for which the 
phase diagram is available. A phase boundary for small 
world to regular transition has also been obtained. The 
network may have small world behaviour even when the 
degree distribution is exponential or stretched exponen- 



tial. Along the a — line, the network retains the small 
world behaviour, consistent with the results of [9] , where 
it was found that D assumed a finite value {^ InN) for 
networks of different sizes for all values of (3. 

It is worth mentioning here that the limiting forms of 
the degree distribution, at extreme values of a, are delta 
functions in nature, but we have restricted our analysis 
to finite values of a. Also, we find that the phase dia- 
gram shows varied features for values of a < for which 
the model corresponds to realistic networks like citation, 
collaboration or social networks. 
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